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Abstract  A  general  process  for  optimization  of  a  sandwich 
panel  to  minimize  the  effects  of  air  blast  loading  is  presented 
here.  The  panel  geometry  consists  of  two  metal  face  plates 
with  a  crushable  honeycomb  or  other  type  of  core.  Opti¬ 
mization  is  necessary  as  there  is  strong  coupling  between 
the  several  variables  and  the  physics,  which  makes  para¬ 
metric  studies  relatively  ineffective.  Virtual  testing  is  used 
to  develop  a  homogenized  model  for  the  stress-strain  curve 
of  the  honeycomb  core,  which  can  be  readily  applied  to 
other  types  of  cellular  core.  The  homogenized  model  has 
been  validated  by  comparison  to  existing  results  as  well 
as  to  results  from  detailed  finite  element  (FE)  models.  A 
design  of  experiments  (DOE)  based  response  surface  opti¬ 
mization  method  in  combination  with  LS-DYNA  is  used  to 
minimize  dynamic  deflection  or  acceleration  of  the  back 
face  plate.  Constraints  on  total  mass  and  on  plastic  strain 
in  the  face  plates  are  imposed.  The  mechanism  of  lower¬ 
ing  the  backface  deflection  is  by  increasing  front  face  plate 
thickness  which  effectively  distributes  the  blast  load  to  a 
larger  area  of  the  core  and  avoids  local  concave  deforma¬ 
tion  of  the  front  face  plate.  Further,  core  depth  is  increased 
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which  increases  panel  stiffness.  For  acceleration  minimiza¬ 
tion,  results  again  produce  a  stiffer  front  face  plate,  but 
accompanied  by  a  sufficiently  soft  core.  The  mechanism  of 
lowering  the  backface  acceleration  is  by  absorbing  energy 
with  low  transmitted  stress.  A  clear  cut  comparison  between 
monolithic  metal  plates  and  sandwich  plates,  for  the  same 
loading  and  failure  criteria,  is  presented  here. 

Keywords  Honeycomb  •  Homogenization  •  Blast  • 
Optimization  •  Sandwich  plates  •  Virtual  testing 

1  Introduction 

Sandwich  panels,  such  as  a  honeycomb  core  with  two  metal 
facing  plates,  are  finding  increasing  use  over  monolithic  or 
solid  plates  in  structural  design  to  withstand  intense  short 
duration  pressure  loads.  Applications  include  protection  of 
land  vehicles,  ships  or  other  structures.  The  cellular  core 
has  the  ability  to  absorb  the  impact  energy  of  the  pres¬ 
sure  pulse  by  undergoing  large  plastic  deformation  at  almost 
constant  nominal  stress.  This  characteristic  of  the  cellular 
core  results  in  significant  reduction  in  the  backface  acceler¬ 
ation  and  hence  mitigates  the  damage  causing  potential  of 
the  blast  impulse.  Though  metal  sandwich  panels  have  been 
used  for  a  long  time  in  aircraft  and  other  light  weight  struc¬ 
tures  to  maximize  the  bending  stiffness  per  unit  density, 
only  recently  have  researchers  investigated  the  possible  use 
of  sandwich  panels  for  blast  protection.  In  sandwich  panels, 
studies  involving  one  or  at  most  two  parameters  at-a-time 
have  been  carried  out.  However,  optimization  with  several 
variables  is  necessary  to  capture  interacting  physics  which 
include:  (1)  too  thin  a  front  face  plate  will  result  in  sig¬ 
nificant  concave  deformation  under  load  that  will  increase 
the  momentum  on  the  structure  which  is  detrimental,  (2) 
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front  and  back  face  plate  thicknesses  must  ensure  proper 
load  transfer  to  the  core  to  enable  crushing,  (3)  the  core 
itself  must  be  stiff  enough  to  minimize  peak  backface  dis¬ 
placement  and  soft  enough  to  crush  and  absorb  energy  to 
minimize  acceleration,  (4)  total  mass  and  space  or  envelope 
constraints  must  be  satisfied,  and  (5)  face  plates  must  be 
thick  enough  to  maintain  integrity. 

In  this  paper,  a  process  is  presented  which  enables  simul¬ 
taneous  consideration  of  several  variables  relating  to  core 
and  face  plate  geometry.  A  key  step  in  the  process  is  a 
validated  approach  to  generating  a  homogenized  model  for 
the  honeycomb  core.  A  design  of  experiments  (DOE)  based 
response  surface  optimization  method  in  combination  with 
LS-DYNA  is  used  to  minimize  dynamic  deflection  or  accel¬ 
eration  of  the  back  face  plate.  Constraints  on  mass  and 
on  plastic  strain  in  the  face  plates  are  imposed.  Further,  a 
clear  cut  comparison  between  monolithic  metal  plates  and 
sandwich  plates  for  the  same  loading  and  failure  criteria  is 
presented  here. 

The  shape  of  a  monolithic  aluminum  plate  to  reduce 
dynamic  displacement  under  blast  was  optimized  in  Argod 
et  al.  (2010)  and  Belegundu  et  al.  (2008).  For  this,  a 
differential  evolution  optimizer  was  coupled  to  FS-DYNA 
to  minimize  the  plate’s  peak  displacement  subject  to  limits 
on  mass  and  on  peak  plastic  strain;  the  shaped  plate  per¬ 
formed  very  well  compared  to  a  flat  plate  of  equal  mass. 
Xue  and  Hutchinson  (2004)  compared  the  performance  of 
sandwich  panels  (such  as  pyramidal  truss  core,  square  hon¬ 
eycomb  and  folded  plate)  to  a  monolithic  plate  of  equal 
weight  for  blast  resistance.  They  found  that  square  hon¬ 
eycomb  and  folded  plate  outperformed  the  pyramidal  truss 
core,  but  all  three  sandwich  panels  were  capable  of  offering 
higher  blast  resistance  compared  to  the  monolithic  plate. 
Their  optimization  study  did  not  consider  material  failure, 
shape  and  certain  core  parameters.  Fleck  and  Deshpande 
(2004)  developed  an  analytical  methodology  to  analyze  the 
dynamic  response  of  metallic  sandwich  beams  subject  to 
both  air  and  water  blasts.  Their  finding  on  the  basis  of  sim¬ 
ple  analytical  formulas  matched  well  with  the  result  from 
Xue  and  Hutchinson’s  (2004)  three-dimensional  FE  calcu¬ 
lations.  Yen  et  al.  (2005)  carried  out  both  experimental  and 
computational  analyses  to  study  the  effect  of  honeycomb 
crush  strength  on  the  dynamic  response  of  a  honeycomb 
core  sandwich  pendulum  system.  The  result  indicated  that 
total  impulse  of  the  system  increased  due  to  concave  defor¬ 
mation  of  the  front  face  plate,  and  that  significant  reduction 
in  maximum  stress  amplitude  propagating  within  the  core 
can  be  achieved  by  suitable  selection  of  honeycomb  mate¬ 
rial  with  proper  crush  strength.  Further,  suitable  shape  of 
the  front  face  plate  can  reduce  the  concave  deformation 
and  hence  the  total  blast  impulse.  Numerical  analyses  car¬ 
ried  out  in  FS-DYNA  using  ConWep  air  blast  function 
validated  the  experimental  result.  Hanssen  et  al.  (2002) 


performed  similar  tests  on  aluminum  foam  core  sandwich 
panels  with  similar  results.  Main  and  Gazonas  (2008)  inves¬ 
tigated  the  uniaxial  crushing  of  a  cellular  sandwich  plate 
subjected  to  air  blast.  The  study  showed  that  the  shock  trans¬ 
mission  can  be  reduced  by  suitably  distributing  the  mass 
among  the  face  plates  and  the  core  for  a  given  mass  of 
the  sandwich.  The  study  considered  the  effect  of  face  plate 
thicknesses  and  core  depth  but  did  not  consider  parame¬ 
ters  related  to  the  core  such  as  core  density.  They  showed 
that  prior  to  densification,  the  core  provides  structural  sup¬ 
port  to  the  front  plate  and  regulates  the  stress  transferred 
to  the  back  face  plate.  Once  the  onset  of  core  densification 
starts,  higher  stresses  are  transferred  to  the  back  face  plate 
(Main  and  Gazonas  2008;  Chi  et  al.  2010).  Zhu  et  al.  (2008, 
2009)  presented  a  limited  optimal  design  study  of  the  hon¬ 
eycomb  sandwich  panel  and  showed  that  there  exists  an 
optimum  core  density  and  core  depth  to  minimize  the  sand¬ 
wich  deflection.  Karagiozova  et  al.  (2009)  states  that  the 
optimum  sandwich  configuration  depends  upon  the  applied 
blast  load,  and  an  optimum  structure  compromises  between 
energy  absorption  of  the  core  and  the  load  transfer  to  the 
back  face  plate  of  sandwich.  For  the  same  load,  the  acceler¬ 
ation  of  the  back  face  plate  also  depends  upon  the  mass  of 
the  back  face  plate. 

It  is  relevant  to  also  cite  papers  that  address  the  stress- 
strain  response  of  honeycomb  sandwich  panels,  as  these 
have  helped  to  validate  portions  of  the  homogenization  pro¬ 
cess  in  this  paper.  Yamashita  and  Gotoh  (2005)  studied 
the  impact  behavior  of  honeycomb  cells  through  numerical 
simulations  and  experiments.  Numerical  simulation  using 
a  single  ‘Y’  cross-sectional  unit  cell  model  predicted  the 
crush  behavior  quite  well  compared  to  experiments  with 
drop  hammer  velocity  of  10  m/s.  Highest  crush  strength 
per  unit  mass  was  obtained  when  cell  shape  is  a  regular 
hexagon.  Wierzbicki  (1983)  derived  a  simple  formula  from 
the  basic  principles  of  material  continuity  and  plasticity  for 
calculating  the  mean  crush  strength  of  metal  honeycombs 
in  terms  of  the  cell  diameter,  foil  thickness  and  the  flow 
stress.  The  derivation  is  given  for  a  general  shape,  and  is 
then  specified  for  a  regular  hexagon  cell.  The  result  from 
this  analytical  solution  is  well  matched  to  the  experimen¬ 
tal  results.  Zhang  and  Ashby  (1992)  analyzed  the  collapse 
behavior  of  the  honeycomb  under  both  axial  compression 
and  in-plane  shear  load.  Buckling,  debonding  and  fracture 
are  identified  as  possible  collapse  mechanisms.  For  flexible 
honeycombs  such  as  those  made  from  Nomex,  buckling  and 
fracture  are  dominant  mode  of  failure  in  simple  axial  com¬ 
pression  test,  but  for  rigid-plastic  honeycombs  (made  from 
aluminum),  buckling  and  plastic  yielding  dominates.  Depth 
of  the  honeycomb  has  no  effect  and  cell  angle  has  little 
effect  on  out-of-plane  strengths  (compressive  and  shear). 
These  strengths  are  highly  sensitive  to  the  density  of  the 
honeycomb.  It  is  also  found  that  out-of-plane  loading  has 
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little  effect  on  in-plane  failure  and  vice  versa.  Wu  and  Jiang 
(1997)  performed  both  quasi-static  and  high  speed  impact 
(up  to  28.14  m/s)  crush  test  on  six  types  honeycomb  cel¬ 
lular  structure.  They  mentioned  that  small  cell  size,  short 
height  honeycomb  made  from  high  strength  material  has 
high  energy  absorbing  capacity. 


2  Homogenization  of  honeycomb  structure 
via  nonlinear  virtual  testing 

It  is  necessary  to  homogenize  the  honeycomb  core,  as  a 
detailed  finite  element  model  of  a  sandwich  panel  with  a 
core  will  require  a  high  density  mesh  to  capture  the  cyclic 
plastic  buckling  or  ‘folding’  deformation  of  the  core  accu¬ 
rately,  entailing  enormous  computing  time  for  each  analysis. 
Moreover,  optimization  involves  iterative  analysis.  Homog¬ 
enization  will  allow  the  honeycomb  to  be  replaced  by  a  3-D 
continuum  structure  which  can  then  be  modeled  using,  say, 
eight-noded  hexahedral  elements.  A  process  is  given  here 
which  can  be  used  to  homogenize  any  structural  concept  for 
the  core.  The  main  steps  are:  (1)  a  unit  cell  is  identified 


and  a  FE  model  of  the  cell  is  developed,  (2)  geometrical 
parameters  associated  with  the  unit  cell  are  identified,  (3) 
nonlinear  virtual  testing  is  carried  out  for  different  values  of 
the  geometrical  parameters,  followed  by  curve  fitting  which 
parametrizes  the  stress-strain  curve  in  terms  of  the  geomet¬ 
rical  parameters,  (4)  validation  of  the  homogenized  model 
both  by  detailed  FE  modeling  of  the  unit  cell  and/or  of  a  test 
specimen  and  by  comparison  to  any  existing  results  in  the 
literature.  Since  the  actual  blast  panel  undergoes  compres¬ 
sion  as  well  as  bending,  validation  must  include  different 
types  of  virtual  tests.  Following  this,  optimization  of  the 
geometry  of  both  the  honeycomb  core  and  of  the  face  plates 
is  carried  out.  These  four  steps  are  described  below. 

2.1  Unit  cell 

Figure  1  shows  the  hexagonal  cell  structure,  its  unit  cell,  and 
its  FE  model.  The  simplest  repeating  unit  in  this  structure  is 
a  ‘Y’  shape,  which  is  taken  to  be  the  unit  cell  (Yamashita 
and  Gotoh  2005).  The  regular  hexagon  cell,  branch  angle  = 
120°,  is  considered  here  as  it  gives  the  highest  crush  strength 
per  unit  mass  (Yamashita  and  Gotoh  2005).  The  unit  cell 


(a)  (b)  (c) 


Fig.  1  Numerical  model  of  honeycomb  core  (a)  Honeycomb  cell  geometry,  (b)  Unit  cell,  (c)  Boundary  conditions  on  the  unit  cell,  (d)  FEM  in 
LSDYNA  and  (e)  Adhesive  in  model 
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Table  1  Material  properties  of 
the  honeycomb  unit  cell  model 

Material 

Density 

(kg/m3) 

Young’s  modulus 
(GPa) 

Yield  Stress 

(MPa) 

Tangent  modulus 
(MPa) 

Poison’s 

ratio 

Foil-AL5052 

2680 

72 

300 

50 

0.34 

Adhesive 

2000 

5 

30 

0 

0.3 

Drop  Hammer-Steel 

288  x  105 

200 

- 

- 

0.24 

has  one  double  wall  and  two  single  walls.  The  double  wall 
consists  of  two  layers  of  foil  glued  together  by  an  adhesive. 
LS-DYNA  is  used  to  carry  out  the  virtual  simulation.  The 
foil  is  modeled  by  quadrilateral  Belytschko-Tsay  shell  ele¬ 
ments,  and  the  0.01  mm  thick  layer  of  adhesive  at  the  double 
wall  is  modeled  by  solid  elements.  Symmetric  boundary 
conditions  are  applied  along  all  the  edges  of  the  foil,  bottom 
areas  are  fixed  and  displacement  load  (crushing)  is  applied 
to  an  external  rigid  surface  via  a  drop  hammer  which  hits 
the  top  areas  and  moves  with  them.  The  main  role  of  the 
top  and  bottom  face  plate  is  to  contain  the  crushed  hon¬ 
eycomb  foil.  In  the  actual  mechanical  test,  a  heavy  steel 


hammer  is  used  to  crush  the  honeycomb.  Material  proper¬ 
ties  are  given  in  Table  1.  To  replicate  the  actual  test,  the 
rigid  surface  is  modeled  using  the  rigid  shell  element  and 
the  mechanical  properties  are  defined  as  that  of  steel  but 
with  a  high  fictitious  density.  AL5052  aluminum  alloy  with 
bilinear  isotropic-hardening  elastoplastic  material  model  is 
used  for  the  foil.  Since  the  yield  and  ultimate  strength  of 
the  AL5052  foil  are  very  close,  bilinear  elastoplastic  mate¬ 
rial  model  with  very  low  tangent  modulus  is  a  reasonable 
approximation.  The  adhesive  is  modeled  as  perfectly  plas¬ 
tic.  Automatic  single  surface  contact  is  applied  to  the  model 
with  sliding  and  sticking  frictional  coefficients  equal  to  0.2 


(a)  (b)  (c)  (d) 


Densification 


0.0  . 

0.0  0.2  0.4  0.6  0.8  1.0 


Volumetric  Strain 

(e) 

Fig.  2  Honeycomb  core  crushing  (a)-(d)  Different  stages  of  honeycomb  unit  cell  crushing  and  (e)  obtained  load  curve  and  its  different  parameters 
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Fig.  3  Variation  of  crush  strength  with  t/D  obtained  from  virtual 
testing 
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Fig.  4  Validation  of  the  crush  strength  obtained  from  the  virtual  test 


and  0.3,  respectively.  These  mechanical  properties  of  the 
foil  and  the  adhesive,  and  friction  coefficient  values  are 
obtained  from  the  literature  (Yamashita  and  Gotoh  2005). 
All  the  tests  are  carried  out  at  a  hammer  speed  of  80  mm/s 
along  the  depth  of  the  honeycomb,  which  can  be  consid¬ 
ered  as  quasi- static  tests  since  impact  tests  involve  very  high 
hammer  speeds  in  the  order  of  m/s. 


2.2  Geometrical  parameters  related  to  unit  cell 

The  following  parameters  of  the  honeycomb  cell  geometry 
were  considered:  foil  thickness  t,  cell  size  D  which  equals 
the  distance  between  opposite  walls  of  the  honeycomb  cell, 
and  core  depth,  h.  Preliminary  FE  runs  show  that  the  stress- 
strain  or  load  curve  does  not  show  any  visible  change  with 
core  depth  h  (Zhang  and  Ashby  1992),  provided  that  it  is 
not  too  small  so  as  not  to  allow  folds  to  occur  during  crush¬ 
ing.  Further  the  ratio  t/D  is  the  defining  parameter  that 
characterizes  the  load  curve.  This  has  also  been  reported 
in  Wierzbicki  (1983).  The  density  is  linearly  proportional  to 
t/D,  as  will  be  given  below. 


2.3  Nonlinear  virtual  testing  and  parametrization  of  the 
stress-strain  curve 

Figure  2  shows  a  typical  crushing  phenomenon  and  load 
curve  obtained  from  the  FE  test.  As  the  hammer  travels, 
buckling  (Yamashita  and  Gotoh  2005;  Zhang  and  Ashby 
1992)  of  the  foil  starts  from  near  the  impact  edge  and  propa¬ 
gates  downward.  Compressive  stress  is  obtained  by  dividing 
the  reaction  force  experienced  by  the  hammer  by  the  trian¬ 
gular  unit  cell  area,  and  volumetric  strain  is  calculated  from 
change  in  core  depth  divided  by  its  original  value.  The  core 


Fig.  5  Detailed  model  of  test  specimen 
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resists  buckling  until  the  peak  stress  point.  Onset  of  buck¬ 
ling  causes  a  sudden  drop  in  the  compressive  stress  (Fig.  2a). 
Compressive  stress  drops  until  the  first  folding  of  the  cell 
wall  is  complete  (Fig.  2b)  and  then  stress  increases.  This 
goes  on,  although  the  peaks  are  very  small  in  comparison 
to  the  first  peak,  until  the  whole  core  is  folded.  The  crush 
strength  (i.e.  the  plateau  stress  level)  is  the  average  of  the 
oscillatory  stress  during  the  cyclic  collapse  of  the  foil.  Once 
the  entire  core  is  folded,  densification  starts  resulting  in  very 
high  compressive  stress.  Large  amount  of  energy  gets  dis¬ 
sipated  through  the  plastic  deformation  of  the  cell  wall  at 
each  folding.  The  energy  absorbed  per  unit  area  of  the  core 
is  essentially  equal  to  the  area  under  the  load  curve  multi¬ 
plied  by  core  depth.  Although  sufficient  care  has  been  taken 
in  approximating  the  load  curve,  it  is  not  possible  to  define 
the  crush  start  and  end  strain  very  accurately. 

By  repeated  virtual  tests  for  different  values  of  t/D ,  and 
using  curve  fitting  (Fig.  3  shows  this  for  crush  strength),  we 
obtain  the  following  equations  that  parameterize  the  load 
curve  in  terms  of  honeycomb  geometry.  Below,  p/  =  den¬ 
sity  of  foil  material  in  kg/m3  and  Sy  =  yield  stress  of  the  foil 
material  (aluminum,  here).  Mass  density  refers  to  the  mass 
per  unit  cell  area  and  per  unit  depth.  Yield  stress  and  density 
refer  to  that  of  AL5052  material.  Crush  start  and  end  strain 
do  not  change  much  with  t/D  and  corresponding  strain  val¬ 
ues  are  taken  as  0.009  and  0.744  respectively.  Final  strain  is 
fixed  at  0.85. 

Mean  crush  stress  (MPa)  =  2296  ( t/D)lA 9 
Peak  strain  =  0.0626  (t/D)  +  0.0035 
Peak  stress  (MPa)  =  845.8  (t/D)  -  0.959 
Young’s  modulus  =  Peak  stress /Peak  strain 
Crush  start  strain  =  0.009 
Crush  end  strain  =  0.744 
Densification  strain  =  0.85 
Densification  stress  (MPa)  =  35592  ( t/D )1,69 

Tensile  stress  cutoff  (MPa)  =  (t/D) 

Mass  density  (kg/m3^  =  (t/D)  (1) 

Relations  in  (1)  provide  a  homogenized  model  for  the 
honeycomb  core. 

2.4  Validation  of  the  homogenized  model 

Several  validation  techniques  are  discussed  here.  Firstly, 
Wierzbicki  (1983)  showed  that  the  mean  crush  strength  may 
be  given  by 

Mean  crush  strength  =  16.565 j(Vd)  ^  (2) 


Fig.  6  Comparison  of  load  curves  from  unit  cell  and  detailed  model 
on  a  virtual  test  specimen.  Numbers  in  the  parenthesis  is  the  average 
crush  strength 


The  analytically  derived  formula  in  (2)  shows  good  agree¬ 
ment  with  the  corresponding  equation  obtained  from  virtual 
testing  in  (1)  as  shown  in  Fig.  4.  Secondly,  the  load  curve 
obtained  from  the  model  in  Fig.  1  is  compared  to  that 
obtained  from  a  detailed  model  of  a  test  specimen  (Fig.  5). 
The  load  curves  match  well  (Fig.  6).  Here,  various  bound¬ 
ary  conditions  have  been  virtually  tested  including  random 
imperfections  in  geometry,  all  producing  good  agreement 
with  the  baseline  test  reported  above. 

Thirdly,  since  the  homogenized  model  is  based  on  crush¬ 
ing,  and  the  blast  panel  is  subject  to  both  bending  and 
crushing,  validation  via  both  bending  and  crushing  virtual 
were  performed.  A  3-point  bend  test  comparing  the  homog¬ 
enized  model  above  with  a  detailed  virtual  test  specimen 
was  carried  out  (Fig.  7).  Provisions  of  ASTM  standards  D 
7250  and  C  393  are  used  (ASTM  2006a,  b).  For  this  test, 
t/D  =  0.00754  is  used.  The  response  matches  well  for 
this  bending  application  as  well  (Fig.  8).  Comparison  of  a 
beam  with  span  of  beam  perpendicular  to  the  ribbon  was 
also  carried  out  with  good  agreement.  Further,  energy  cal¬ 
culations  were  also  in  good  agreement.  Many  details  related 


Fig.  7  Honeycomb  ribbon  orientation  parallel  to  span  of  sandwich 
beam-detailed  model 
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Fig.  8  Maximum  deflection  of 
the  sandwich  vs.  contact  load 
when  honeycomb  ribbon  is 
parallel  to  the  span  of  sandwich 
beam 


to  loading,  boundary  conditions  and  material  properties  in 
this  beam  test  are  omitted  for  brevity  (Singh  2011),  as  the 
focus  is  on  showing  a  process  by  which  we  can  be  confident 
of  using  the  homogenized  model  for  optimization. 


3  Optimization  problem  formulation 

Optimal  design  of  a  square  honeycomb  core  sandwich  panel 
subjected  to  air  blast  loading  is  considered  (Fig.  9).  The 
panel  is  freely  suspended.  The  role  of  the  stiffener  at  the  top 
is  to  impose  high  inertia  to  the  back  face  plate  and  hence 
the  sandwich,  permitting  compression  of  the  core.  This  set 
of  boundary  conditions  are  motivated  by  design  problems  in 
vehicle  protection  and  approximate  commonly  used  exper¬ 
imental  fixtures.  A  study  of  different  boundary  conditions 
on  monolithic  plates,  in  contrast  with  sandwich  panels  as 
considered  here,  was  previously  conducted  in  Argod  et  al. 
(2010). 

In  general,  design  variables  include  thicknesses  of  face 
plates,  core  depth,  core  t/D  or  equivalently  core  density, 
and  bulge  magnitudes  of  the  face  plates  (Fig.  10).  Thus,  x  = 

{ tb ,  tf ,  h ,  t/D ,  Sb,  Sf}. 


Two  separate  objectives  are  considered  for  minimization: 
one  is  the  peak  displacement  magnitude  and  the  other  is  the 
rigid  body  acceleration,  of  the  back  face  plate.  Constraints 
are:  limits  on  total  mass  and  on  plastic  strain  limit  in  the 
face  plates,  the  latter  to  ensure  structural  integrity.  Thus,  we 
have  the  problem 


Minimize 

max  8b 

(case  —  1) 

ma  xotb 

(case  —  2) 

subject  to 

£pj  ^  Spmax 

for  each  element  j  in  the 
face  plates 

M  <  MmSLX 

mass  limit 

xL  <  x  <  xu 

bounds  on  variables 

In  case-1,  since  the  sandwich  model  is  not  constrained,  the 
maximum  Z-deflection  of  the  back  face  plate  is  obtained 
by  subtracting  the  rigid  body  displacement  of  the  stiffener 
as  8b  =  <$z_ max  —  ^stiffener-  Displacements  along  the  x- 
and  y-  direction  are  not  significant  and  are  not  considered. 
The  displacement  is  a  function  of  time,  and  the  peak  value 
is  monitored.  In  case-2,  the  absolute  maximum  rigid  body 
acceleration  of  back  face  plate  along  z-direction  is  consid¬ 
ered.  Like  displacement,  accelerations  along  the  x-  and  y- 
direction  are  not  significant.  Rigid  body  acceleration  refers 


Fig.  9  (a)  Honeycomb 
sandwich  panel  (b)  Exploded 
view  of  the  honeycomb  core 
sandwich  model  used  for 
optimization  study 
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Fig.  10  Optimization  parameters  for  the  sandwich  optimization 


to  net  force  divided  by  net  mass.  The  plastic  strains  in  the 
face  plates  increase  with  simulation  time  until  a  plateau  is 
reached  and  this  saturated  or  maximum  value  is  considered. 
The  maximum  value  of  the  all  elements  in  the  face  plates  is 
monitored.  The  mass  M  is  the  sum  of  the  mass  of  the  front 
and  back  face  plates  and  of  the  core  as  M  =  Mf  +  Mb  +  Mc . 
Mass  of  the  stiffener  remains  constant  during  optimization 
and  is  not  included  in  M. 

Optimization  is  carried  out  in  the  following  sequence:  (I) 
optimization  for  case-1  considering  only  {tb,  tf,  h ,  t/D} 
as  variables,  (II)  optimization  for  case-1  considering  all  six 
variables  {tb,  tf,  h ,  t/D,  Sb ,  Sf},  and  (III)  optimization  for 
case-2  considering  only  { tb,tf,h ,  t/D}. 

For  brevity  in  the  figure  captions,  we  refer  to  {tb,  tf, 
h ,  t/D} -optimized  panel  as  ‘size-optimized  panel’,  and  {tb, 
tf ,  h ,  t/D ,  sb ,  Sf]  -optimized  panel  as,  ‘size+shape  opti¬ 
mized  panel’ .The  variable  t/D  is  a  parameter  that  defines 
the  stress-strain  curve  which  is  input  into  LS-DYNA.  The 
remaining  five  variables  {tb,  tf,  h ,  Sb,  Sf}  require  a 
change  in  the  coordinates  of  the  nodes  in  the  3D  finite 
element  model.  That  is,  they  affect  the  shape  of  the  struc¬ 
ture.  The  key  equation  to  implement  shape  optimization  is 
(Belegundu  and  Rajan  1988) 

Ndv 

G  (x)  =  Goriginal  +  ^  (4) 

i= 1 

where  G  is  a  grid  point  coordinates  vector,  representing  x-, 
y-,  z-  coordinates  of  all  nodes  in  the  model.  Each  repre¬ 
sents  the  amplitude  of  a  ‘permissible  shape  change  vector’ 
or  what  is  commonly  called  a  ‘velocity  field’  vector  q/v\ 
Velocity  fields  have  nothing  to  do  with  actual  velocities  of 
the  model  under  loading.  Vectors  { }  are  generated  just 
once  in  the  optimization  procedure.  Goriginal  is  the  current 
(flat)  shape.  Visualization  of  a  { qz }  is  identical  to  visual¬ 
ization  of  a  displacement  field  in  finite  elements:  { q^ }  is 
multiplied  by  a  magnification  scalar  and  added  to  the  cur¬ 
rent  grid  to  obtain  a  displaced  grid,  except  that  here  the 
displaced  grid  represents  a  new  shape  and  is  called  a  basis 
shape. 

Velocity  fields  associated  with  variables  tb,  tf ,  h  are 
straightforward:  nodes  are  moved  in  the  +/—  z-direction, 


maintaining  equal  spacing,  to  result  in  the  desired  thick¬ 
nesses.  For  variables  Sb  and  Sf ,  which  are  the  amplitudes 
of  basis  shapes  that  correspond  to  bulges  in  the  face  plates 
as  shown  in  Fig.  10,  a  procedure  is  needed  to  generate  q* . 
The  velocity  fields  q*  are  generated  here  by  first  applying 
a  dummy  load  on  each  of  the  plate  surfaces  and  using  ana¬ 
lytical  expressions  in  Timoshenko  and  Woinowsky-Krieger 
(1959),  Eq.  145,  page  142,  to  obtain  the  corresponding 
nodal  displacements  on  all  nodes  on  the  surface.  Then,  equal 
spacing  in  the  z-direction  is  used  to  complete  the  definition 
of  q*.  A  square  portion  (1.016  m  x  1.016  m)  at  the  cen¬ 
ter  of  the  face  plates  is  taken  as  the  domain  for  applying 
velocity  fields.  Asymmetric  shape  variations  are  not  nec¬ 
essary  in  this  problem  owing  to  centrally  located  charge. 
Eower  bounds  Sb  and  s /  are  zero  implying  that  only  con¬ 
vex  (outward)  bulges  are  allowed  for  shape  changes,  as  this 
has  found  to  be  beneficial  in  deflecting  the  waves  in  the 
monolithic  plate  (Argod  et  al.  2010;  Belegundu  et  al.  2008). 

As  variables  x  =  {tb,  tf,  h,  t/D ,  Sb,  Sf}  are  changed 
for  each  sampling  point,  the  grid  point  coordinates  G,  the 
load  curve  from  (1),  and  thickness  values  are  updated,  an 
ES-DYNA  input  file  is  then  written,  and  an  analysis  is  car¬ 
ried  out  to  evaluate  the  various  functions  in  the  optimization 
problem.  The  process  is  schematically  given  as: 


Use  current  values  of  [tb,  tf,  h,  Sb,  Sf }  at  sampling  point 
xk,  construct  G  using  { qz } ;  from  current  value  of  {t/D}, 
construct  the  stress-strain  curve  data  using  (1)  — ►  Write 
input  file  and  perform  LS-DYNA  analysis  — ►  Evaluate 
objective  and  constraint  functions  — ►  Create  response 
equations  — ►  Run  optimizer 


The  DOE  response  surface  method  is  used  to  solve  this 
problem,  implemented  using  Design  Expert,  a  commer¬ 
cially  available  software.  Using  the  response  equations, 
optimization  is  carried  out  by  fmincon,  a  gradient  based 
optimizer  in  MATLAB  optimization  toolbox.  Central  com¬ 
posite  face  centered  (CCF)  method  is  adopted  to  create 
design  points.  It  uses  three  levels  for  each  factor.  Lower 
design  limits  are  decided  in  such  a  way  that  FE  evaluation 
is  feasible  at  all  sampling  points  and  to  avoid  high  aspect 
ratios  in  the  hexahedral  elements.  Also,  appropriate  upper 


Table  2  Design  limits  for  {tb,  tf,  h,  t/D} -optimization,  min.  8b 


Design  limits  Lower  limit  Upper  limit 


Front  face  plate  thickness  (mm),  tf 

4.4 

18 

Core  depth  (mm),  h 

280 

500 

Foil  thickness/cell  size,  t/D 

0.018 

0.046 

Back  face  plate  thickness  (mm),  tb 

4.4 

10 

4^  Springer 
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Table  3  Optimization  results  for  {tb,  tf,  h,  t/D] -optimization  with 
varying  mass  limits,  min.  8b  (units:  mm) 


Sandwich  Optimized  parameters  £Pmax 


mass  (kg) 

tf 

h 

tb 

(t/D) 

5b  (Obj) 

130 

8.07 

300.6 

4.4 

0.0252 

24.3 

0.0384 

140 

7.90 

359.9 

4.4 

0.0238 

19.96 

0.0373 

150 

7.55 

401.6 

4.4 

0.0240 

16.92 

0.033 

160 

8.4 

390.2 

4.4 

0.0263 

14.65 

0.023 

170 

9.12 

388 

4.4 

0.0282 

12.43 

0.017 

design  limits  are  taken  so  that  better  response  fit  is  possi¬ 
ble  and  at  the  same  time  giving  more  design  space  to  the 
optimizer  to  find  the  minima.  Design  limits  for  the  different 
cases  are  given  in  Tables  2,  3  and  4  below.  Some  results 
have  been  cross-checked  by  using  the  differential  evolution 
(DE)  technique  which  is  computationally  more  expensive. 


4  FE  modeling  using  LS-DYNA 

In  this  study,  aluminum  AL5052  is  considered  for  the  entire 
sandwich  such  that  results  can  be  easily  compared  with 
AL5052  monolithic  plate  considered  earlier  (Argod  et  al. 
2010;  Belegundu  et  al.  2008).  In  the  monolithic  plate  prob¬ 
lem,  various  boundary  conditions  were  considered  giving 
similar  optimum  shapes  (Argod  et  al.  2010).  The  model  is 
free  to  move  in  space,  which  approximates  to  some  extent 
commonly  used  experimental  fixtures.  As  a  consequence  of 
this,  the  back  face  plate  is  not  restrained  and  can  deform 
freely  without  creating  unrealistically  high  plastic  strain. 
The  role  of  the  stiffener  at  the  top  is  to  impose  high  inertia 
to  the  back  face  plate  and  hence  the  sandwich,  permitting 
compression  of  the  core.  Mass  of  the  stiffener  is  1,850  kg. 
High  fictitious  density  is  defined  for  the  stiffener.  The  con¬ 
tacts  between  the  face  plate  and  the  core,  and  between  the 


Table  4  Design  limits  for  {tb,  tf,  h, 
8b 

t/D ,  Sb,  Sf] -optimization,  min. 

Design  limits 

Lower  limit 

Upper  limit 

Front  face  plate  thickness  (mm),  tf 

4.4 

10 

Core  depth  (mm),  h 

150 

500 

Foil  thickness/cell  size,  t/D 

0.018 

0.046 

Back  face  plate  thickness  (mm),  tb 

4.4 

10 

Front  face  plate  bulge  (mm),  5 / 

10 

100 

Back  face  plate  bulge  (mm),  Sb 

0 

100 

back  face  plate  and  the  stiffener  are  defined  using  CON¬ 
TACT  features  and  *TIED_SURFACE_TO_SURFACE_ID 
card.  After  a  mesh  convergence  study  a  38  x  38  x 
2  and  28  x  28  x  3  element  mesh  is  taken  in  x-y 
plane  for  the  face  plates  and  the  core,  respectively.  The 
*MAT_PLASTIC_KINEMATIC  material  model  is  used  for 
the  face  plates  and  *MAT_CRUSHABLE_FOAM  material 
model  is  used  for  the  core.  Material  properties  used  for  the 
face  plate  are  same  as  that  of  the  foil  used  for  honeycomb 
virtual  testing  (Table  1).  The  mechanical  properties  of  the 
core  are  defined  in  (1),  which  are  functions  of  t/D.  Poisson’s 
ratio  is  taken  as  zero  for  the  core.  The  *LOAD_BLAST 
input  parameters  used  are  equivalent  TNT  mass  (8  kg),  type 
of  blast  (air  blast- spherical  charge),  load  curve,  charge  loca¬ 
tion  (0,  0,  —0.4064  m).  The  exposed  bottom  surface  of  the 
plate  is  defined  as  the  area  on  which  blast  load  applied.  A 
low  value  of  damping  ~  le-9  and  a  reduced  tssfac  help  in 
smoother  convergence  in  LS-DYNA. 

5  Optimization  results 

Results  are  discussed  here  for  different  mass  limits  for  the 
sandwich,  while  detailed  response  is  given  for  the  150  kg 
mass  limit.  The  optimization  study  is  carried  out  for  a  fixed 
amount  of  charge,  viz.  8  kg  TNT.  When  the  DOE  response 
surface  optimizer  provides  an  optimum  set  of  design  vari¬ 
ables,  LS-DYNA  is  executed,  and  the  resulting  response 
values  are  used  in  the  tables  below. 

5 . 1  Optimization  results  for  minimum  backface 
displacement,  8b 

We  note  that  8b  refers  to  the  peak  in  the  displacement-time 
response  of  the  backface  relative  to  the  stiffener.  Table  5 
shows  the  {tb,  tf,  h,  t/D} -optimized  sandwich  panel  param¬ 
eters  for  different  mass  limits.  The  £pmax  is  always  at 


Table  5  Optimization  results  for  {tb,  tf,  h,  t/D ,  Sb,  s  /} -optimization 
with  varying  mass  limits,  min.  8b  (units:  mm) 


Sandwich  Optimized  parameters  £Pmax 


mass  (kg) 

tf 

h 

tb 

(t/D) 

Sf 

Sb 

5b 

(Obj) 

130 

4.4 

150 

4.4 

0.0364 

42.41 

0 

20.76 

0.031 

140 

4.4 

150 

4.4 

0.0371 

52.69 

0 

17.04 

0.020 

150 

4.4 

150 

4.4 

0.0377 

62.96 

0 

14.4 

0.014 

160 

4.4 

150 

4.4 

0.0384 

73.24 

0 

12.25 

0.012 

170 

4.4 

150 

4.4 

0.039 

83.52 

0 

10.73 

0.001 
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Table  6  Design  limits  for  {tb,tf,h,  t/D} -optimization,  min.  ab 


Design  limits 

Lower  limit 

Upper  limit 

Front  face  plate  thickness  (mm),  tf 

15 

25 

Core  depth  (mm),  h 

120 

400 

Foil  thickness/cell  size,  t/D 

0.00754 

0.025 

Back  face  plate  thickness  (mm),  tb 

8 

25 

the  center  of  the  front  face  plate.  With  increase  in  mass, 
8b  decreases.  The  optimizer  produces  thicker  and  denser 
(higher  t/D)  core.  Mass  is  always  active  and  the  plastic 
strain  is  active  only  for  130  kg  and  140  kg. 

Table  6  shows  the  {tb,  tf,  h,  t/D ,  Sb,  Sf } -optimized  sand¬ 
wich  panel  parameters  for  different  mass  limits.  The  spmax 
is  always  in  the  back  face  plate.  With  increase  in  mass,  8b 
decreases.  The  optimizer  produces  a  denser  or  higher  t/D 
core.  The  %  is  always  at  lower  limit  (4.4  mm)  for  the  rea¬ 
son  mentioned  in  Section  6.  The  tf  is  always  at  lower  limit 
(4.4  mm)  as  the  front  face  plate  thickness  re-emerges  in  the 
bulge.  The  bulge  also  reduces  the  effective  blast  load  on  the 
panel  by  deflecting  the  blast.  The  bulge  also  provides  higher 
stiffness  to  the  panel,  compensating  for  the  loss  in  stiffness 
due  to  lower  core  depth.  With  increase  in  mass,  t/D  and  s / 
increase  consistently.  Mass  is  always  active  and  the  plastic 
strain  is  always  inactive. 

Now,  focusing  on  the  150  kg  mass  limit,  Fig.  11  shows 
the  performance  of  the  optimized  sandwich  panel  and 
how  it  compares  with  an  all-aluminum  monolithic  plate. 
A  non-optimized  uniformly  thick  aluminum  panel  is  also 
included  in  the  comparison  for  reference  purposes.  The  cor¬ 
responding  backface  deflection-time  responses  are  shown  in 


Fig.  12a,  plastic  strains  in  Fig.  12b  and  Z -momentum  in 
Fig.  12c.  Compared  to  the  uniformly  thick  aluminum  panel, 
the  sandwich  panel  shows  significant  reduction.  Shape 
optimization  of  the  sandwich  adds  to  this  improvement. 
However,  a  shape-optimized  all-aluminum  panel  scored  a 
victory  as  compared  to  the  sandwich  (11.5  mm  to  14.4  mm 
deflection).  A  physical  explanation  for  this  is  given  subse¬ 
quently.  In  the  sandwich,  the  mass  fraction  in  the  front  face 
plate,  core  and  the  back  face  plate  are  0.20,  0.  68  and  0.12, 
respectively,  for  the  {tb,  tf ,  h,  t/D} -optimized  panel  and 
0.48,  0.40  and  0.12,  respectively,  for  the  {tb,  tf ,  h,  t/D ,  Sb, 
Sf } ,  optimized  panel. 


6  Mechanism  causing  the  improvement  in  backface 
deflection 

In  the  {tb,  tf,  h,  t/D}- optimized  panel,  the  optimizer 
has  chosen  a  larger  core  depth  which  increases  the  overall 
stiffness  of  the  sandwich  and  lowers  8b.  The  tb  is  always 
at  lower  limit  (4.4  mm)  as  it  is  sufficient  to  keep  spmax 
below  the  strain  limit  (0.038).  The  tf  is  decided  based  on 
three  things,  viz.  to  keep  spmax  below  the  strain  limit,  to 
provide  sufficient  stiffness  for  the  front  face  plate  to  trans¬ 
fer  the  blast  pressure  load  to  a  larger  area  of  the  core,  and 
to  reduce  local  concave  deformation  of  the  front  face  plate 
which  will  increase  the  impulse. 

Compared  to  the  {tb,tf,h,  t/D} -optimized  panel  above, 
the  {tb,  tf,  h,  t/D ,  sb,  Sf} -optimized  panel  has  further 
reduced  the  deflection  by  14.9  %.  This  may  be  attributed 
to  the  following  mechanisms.  A  lesser  core  depth  has  been 
produced  but  with  a  bulge  in  the  front  face  plate.  The  core 
is  denser.  The  bulge  deflects  the  blast  wave  and  reduces 


Fig.  11  Optimized  sandwich 
panel  and  shape  optimized 
monolithic  plate  of  150  kg  mass 
for  minimizing  8b  (stiffener  not 
shown) 
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Fig.  12  Comparison  for 
minimum  8b  for  150  kg  mass 
limit:  (a)  8b,  (b)  spmax  and  (c) 
total  Z-momentum.  Note:  by 
solid  plate  is  meant  a  monolithic 
all-metal  plate 
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the  impulse  imparted  to  the  structure,  and  it  also  increases 
the  moment  of  inertia  of  the  sandwich  at  its  center  where 
blast  load  is  maximum.  The  bulge  also  makes  the  front  face 


plate  stiffer,  which  helps  to  transfer  the  blast  load  to  a  large 
area  of  the  core  thereby  reducing  local  deformation  at  the 
center. 
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Table  7  Optimization  results  for  {tb,  tf,  h,  t/D] -optimization  with 
varying  mass  limits,  min.  a t  (units:  mm,  m/s2) 


Sandwich 
mass  (kg) 

Optimized  parameters 

Spmax 

tf 

h 

tb 

(t/D) 

ab  (Obj) 

130 

17.94 

194.93 

10.77 

0.00754 

2.77  x  104 

0.0355  (f) 

140 

17.79 

196.84 

13.4 

0.00754 

2.22  x  104 

0.0360  (f) 

150 

17.63 

199.22 

16.02 

0.00754 

1.84  x  104 

0.0371  (f) 

160 

17.45 

202.16 

18.64 

0.00754 

1.57  x  104 

0.0388  (f) 

170 

17.45 

205.2 

21.1 

0.00754 

1.37  x  104 

0.0388  (f) 

In  the  all-aluminum  monolithic  plate,  more  mass  is  avail¬ 
able  to  go  in  terms  of  bulge,  whereas  in  the  sandwich, 
mass  is  distributed  between  face  plates  and  core.  Higher 
the  bulge  height  towards  the  blast,  greater  is  the  reduc¬ 
tion  in  impulse  and  hence  greater  the  reduction  in  back- 
face  deflection.  This  aspect  and  higher  rigidity  makes  the 
monolithic  plate  slightly  better  than  the  sandwich  panel  as 
shown  in  Fig.  11a. 

6. 1  Optimization  results  for  minimum  backface 
acceleration,  a^ 

Table  7  summarizes  the  {4,  tf,  h ,  t/D} -optimized  results 
for  different  mass  of  the  sandwich.  With  increase  in  mass, 
at,  decreases.  It  shows  that  both  h  and  tt,  increases  with 
sandwich  mass,  but  tf  does  not  vary  much.  The  core  density 
(proportional  to  t/D)  always  remains  at  the  lower  bound. 
The  Spmax  is  very  close  to  the  limiting  value  of  0.038  for 
all  the  mass  limits.  It  can  be  observed  that  with  increase  in 
sandwich  mass,  the  optimizer  adds  mass  more  to  the  back 
face  plate  than  to  the  core;  at  not  only  depends  upon  the 
force  transmitted  to  the  back  face  plate  but  also  on  the  back 
face  plate  mass  (Main  and  Gazonas  2008;  Karagiozova  et  al. 
2009). 


Focusing  on  the  150  kg  mass  limit,  Fig.  13  shows  the 
performance  of  the  optimized  sandwich  panel  and  how  it 
compares  with  an  all-aluminum  monolithic  plate.  The  corre¬ 
sponding  acceleration- time  responses  are  shown  in  Fig.  14a. 
{ tb ,  tf,  h,  t/D] -optimized  sandwich  panel  results  in  84  % 
reduction  in  a &  compared  to  a  uniformly  thick  monolithic 
plate  of  equal  mass.  The  location  of  the  spmax  is  at  the  center 
of  the  front  face  plate. 

Now,  comparing  the  optimized  sandwich  to  an  optimized 
all-aluminum  monolithic  panel,  the  sandwich  scores  a  clear 
victory  (compare  1.84  x  104  to  8.88  x  104  m/s2).  Recall 
that  the  victor  was  reversed  for  8t,  minimization.  The  mono¬ 
lithic  plate  has  no  mechanism  for  energy  absorption.  The 
mass  fraction  in  the  front  face  plate,  core  and  the  back  face 
plate  are  0.47,  0.11  and  0.42  respectively. 


7  Mechanism  causing  the  improvement  in  backface 
acceleration 

A  paper  by  Avalle  et  al.  (2001)  explains  the  relation  between 
energy  absorption  and  transmitted  force.  We  note  that  a &  is  a 
direct  measure  of  the  transmitted  force,  and  that  core  density 
p  and  t/D  are  linearly  related.  In  Fig.  15,  the  area  under 
the  three  curves  /  a  d  s  are  the  same.  The  area  represents 
energy  per  unit  volume,  and  at,  denotes  the  corresponding 
stress  value,  which  is  here  the  compressive  stress  in  the  core 
elements  in  contact  with  the  back  face  plate.  For  low  density 
core  pi ,  the  deformation  is  high  and  the  stress-strain  state  is 
in  the  densification  zone  with  high  stress.  For  high  density 
core  P3,  the  energy  is  absorbed  with  low  deformation  and 
high  stress.  The  medium  density  foam  p2  is  optimum  as  at, 
is  lowest.  In  fact,  an  efficiency  index  has  been  defined  in 
(Avalle  et  al.  2001)  as 
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Fig.  13  Optimized  sandwich 
panel  and  shape  optimized 
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Fig.  14  Comparison  for 
minimum  ab  for  150  kg  mass 
limit:  (a)  ab,  (b)  spmax  and  (c) 
total  Z-momentum 
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Higher  efficiency  implies  that  the  core  absorbs  a  given 
amount  of  energy  with  lesser  transmitted  force.  A  look  at 
the  results  show  that  optimization  has  exploited  this  physics: 


t/D  and  core  depth  have  been  adjusted  so  that  even  though 
the  maximum  strain  in  front  face  elements  have  gone  to 
the  densification  region  (strain  >  0.744),  the  strain  in  the 
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Fig.  15  Stress  (at  the  backface)  for  different  density  foams  corre¬ 
sponding  to  same  energy  absorption 


backface  elements  are  well  within  the  crushing  zone  avoid¬ 
ing  complete  densification.  The  t/D  has  been  pushed  to  its 
lower  bound.  Further,  a  comparison  of  r\  between  the  opti¬ 
mized  results  for  a &  vs.  8^  minimization  (Table  8)  clearly 
shows  that  the  optimizer  has  done  its  task:  when  the  objec¬ 
tive  is  ab,  a  more  efficient  core  is  designed.  Of  course, 
optimization  also  adjusts  the  mass  distribution  between  face 
plates  and  core,  subject  to  constraints  on  sandwich  mass, 
plastic  strain  in  plates  and  lower  bounds  on  t/D. 


8  Conclusions 

A  process  for  optimizing  sandwich  panels  to  mitigate  blast 
loading  is  presented.  While  honeycomb  core  has  been  used 
in  the  study,  the  process  is  also  applicable  to  other  types 
of  cellular  core.  Two  independent  design  objectives,  back- 
face  deflection  of  the  plate  and  backface  acceleration  are 
minimized,  subject  to  mass  and  plastic  strain  constraints. 
The  optimization  is  carried  out  using  DOE  response  surface 
methodology.  LS-DYNA  is  used  for  finite  element  simu¬ 
lations.  Virtual  testing  is  used  to  develop  a  homogenized 
model  for  the  stress-strain  curve  of  the  honeycomb  core, 
and  this  model  has  been  validated  by  comparison  to  existing 
results  as  well  as  to  detailed  FE  model  results. 

The  mechanism  of  lowering  the  backface  deflection  is 
by  increasing  front  face  plate  thickness  which  effectively 


distributes  the  blast  load  to  a  larger  area  of  the  core  and 
avoids  increase  in  impulse  stemming  from  local  concave 
deformation  of  the  front  face  plate.  Further,  core  depth 
is  increased  which  increases  panel  stiffness.  In  the  sand¬ 
wich,  the  mass  fraction  in  the  front  face  plate,  core  and 
the  back  face  plate  are  0.20,  0.  68  and  0.12,  respectively, 
for  the  {4,  tf ,  h ,  t/D) -optimized  panel.  Interestingly,  for 
the  same  mass,  the  shape-optimized  monolithic  panel  is 
more  effective  than  the  optimized  honeycomb  core  sand¬ 
wich  panel  (11.5  mm  vs.  14.4  mm).  In  the  all-aluminum 
monolithic  plate,  more  mass  is  available  to  go  in  terms  of 
bulge,  whereas  in  the  sandwich,  mass  is  distributed  between 
face  plates  and  core.  Higher  the  bulge  height  towards  the 
blast,  greater  is  the  reduction  in  impulse  and  hence  greater 
the  reduction  in  backface  deflection.  This  aspect  and  higher 
rigidity  makes  the  monolithic  plate  slightly  better  than  the 
sandwich  panel  as  shown  in  Fig.  11a. 

Considering  acceleration  minimization,  results  produce  a 
stiffer  front  face  plate,  which  helps  to  distribute  the  crush¬ 
ing  load  to  a  wider  region  of  the  core,  and  a  soft  core  by 
reducing  t/D  to  its  minimum  value.  The  mechanism  of  low¬ 
ering  the  backface  acceleration  is  by  absorbing  energy  with 
low  transmitted  stress.  In  this  case,  honeycomb  core  sand¬ 
wich  panel  proves  to  be  significantly  more  effective  than  an 
optimized  monolithic  panel  (1.84e4  m/s2  vs.  8.88e4  m/s2). 
With  increase  in  sandwich  mass,  the  optimizer  adds  mass 
more  to  the  back  face  plate  rather  than  the  core;  accelera¬ 
tion  not  only  depends  upon  the  force  transmitted  to  the  back 
face  plate  but  also  on  the  back  face  plate  mass.  In  the  sand¬ 
wich,  the  mass  fraction  in  the  front  face  plate,  core  and  the 
back  face  plate  are  0.47,  0.  11  and  0.42,  respectively. 

As  noted  earlier,  Yen  et  al.  (2005)  had  carried  out  an 
experimental  and  computational  study,  and  concluded  that 
significant  reduction  in  maximum  stress  amplitude  propa¬ 
gating  within  the  core  can  be  achieved  by  suitable  selection 
of  honeycomb  material  with  proper  crush  strength.  This 
observation  has  been  borne  out  by  the  results  here,  as 
optimization  with  FS-DYNA  achieves  a  proper  balance, 
accounting  for  different  interacting  physics.  A  general  and 
efficient  design  process  has  been  presented. 

Fracture  of  the  face  plates  themselves  has  been  avoided 
by  imposing  a  limit  on  maximum  plastic  strain  during  opti¬ 
mization.  However,  fracture  causing  debonding  between  the 
core  and  the  face  plates  has  not  been  considered  in  this 
study.  Finally,  the  work  here  can  be  extended  to  include 
strain  rate  sensitivity  in  the  homogenization  model. 


Table  8  Comparison  of 
honeycomb  core  energy 
absorption  efficiency 


Case 

t/D 

f  <j  ds  (Pa) 

<7b  (average)  (Pa) 

Y] 

acceleration_  min  (a%  =  1.84  x  104  ms-2) 

0.00754 

702377.6 

1.58  x  106 

44.5  % 

deflection_min  (a%  =  7.25  x  104  ms-2) 

0.024 

569514.2 

5.59  x  106 

10.2  % 
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